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BASIC  PRINCIPLES  OF  ORGANOMETALLIC  VAPOR  PHASE  EPITAXY 

G.B.  Stringfellow 
University  of  Utah 
Salt  Lake  City.  Utah  841 12 

ABSTRACT 

Crystal  grovvth  processes  in  general  and  epitaxy  in  particular  are  often 
discussed  in  terms  three  disciplines:  thermodynamics,  mass  transport  and 
hydrodynamics,  artd  chemical  reaction  kinetics.  This  paper  will  concentrate  on 
two  of  these,  the  thermodynamic  and  kinetic  aspects  of  organometallic  vapor 
phase  epitaxy  (OMVPE).  Three  major  influences  of  thermodynamics  will  be 
discuss^:  1)  Thermodynamics  defines  the  driving  force  arrd  hence  the  upper 
limit  of  growth  rate.  This  occurs  only  when  all  reactants  in  the  system  are 
allowed  to  equilibrate  with  the  substrate.  2)  Thermodynamics  often  controls 
stoichiometry  and  the  solid  composition  of  alloys.  An  understanding  of 
thermodynamic  and  kinetic  con^raints  leads  to  the  ability  to  grow  metastable 
alloys.  3)  The  driving  force  for  phase  separation  and/or  ordering  into  natural 
superlattice  structures  during  growth  is  also  governed  by  therm^ynamics.  The 
actual  ordered  structures  observed  are  dependent  on  the  surface  kinetics.  This 
aspect  of  kinetics  will  be  addressed  in  addition  to  the  kinetics  of  both 
homogeneous  aixi  heterogeneous  chemical  reactions  occurring  for  a  variety  of 
precursor  molecules.  Each  of  these  topics  is  addressed  in  terms  of  fundamental 
concepts,  with  examples  from  recent  research  on  the  OMVPE  growth  of  lll/V 
semiconductor  com(X)unds  and  alloys. 

1.  INTRODUCTION 

Epitaxial  growth  processes,  in  general,  and  OMVPE,  in  particular,  are  highly 
complex  processes.  Indeed,  early  crystal  growth  studies  were  largely  empirical, 
giving  crystal  growth  the  appearance  of  an  art.  This  is  partly  because  of  the 
complex,  multicomponent,  multiphase  systems  that  are  rK>rmally  of  interest  and 
partly  because  the  process  is  dynamic  and  inhomogeneous  phases  are 
inherent.  In  an  effort  to  systematically  study  and  understand  such  a  complex 
system  the  fundamental  processes  occurring  during  epitaxial  growth  are 
commonly  subdivided  into  hydrodynamics  and  mass  transport,  the  kinetics  of 
chemical  reactions  occurring  homogeneously  in  the  gas  p^ase  and 
heterogeneously  at  the  surface,  and  thermodynamics.  We  will  concentrate  on 
thermo^namics  and  kinetics  in  this  paper.  The  hydrodynamic  aspects  of 
OMVPE  will  be  addressed  in  the  following  paper. 


Thermodynamic  aspects  of  epitaxial  growth  are  in  many  ways  the  most 
important.  This  is  especially  true  when  the  grov4h  rate  is  very  slow.  At  low 
growth  rates  and  relatively  high  temperaUires.  the  chemical  reaction  kinetics 
play  less  of  a  role  than  in  very  rapid  crystal  growth  processes.  In  the  limit  of 
infinitely  slow  growth  rates  themiodynamics  defines  the  concentrations  in  the 
vapor  and  solid  phases.  Thermodynamics  also  determines  the  driving  force  for 
any  crystal  grovv^  process,  heixe  defining  the  maximum  growth  rate. 
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Thermodynamics  is  also  used  to  predict  solid  composition  for  many  growth 
conditions.  This  includes  not  only  alloy  composition,  but  also  solid 
stoichiometry,  incorporation  of  impurities,  separation  into  several  solid  phases, 
artd  the  occurrence  of  ordered  superlattice  structures  in  the  solid.  Thus,  the 
thermodynamic  aspects  of  epitaxy  must  be  urxJerstood  before  proceeding  to  the 
kinetic  aspects  of  growth  that  frequently  control  growth  rate  and.  in  many 
situations,  affect  solid  composition,  the  phases  af^earing.  and  the  occurrence 
of  phase  separation  and  ordering.  Increasingly,  an  understanding  of  the  basic 
aspects  of  epitaxy  has  allowed  a  departure  from  the  empirical  approach  to 
crystal  growth. 

2.  THERMODYNAMIC  TREATMENT  OF  OMVPE 

The  equilibrium  state  for  a  two  phase,  a  +B.  system  is  defined  in  terms  of  the 
chemical  potentials. 

Mi®  =Mi®  .  (1) 

where  the  subscript  i  indicates  the  ith  component  and  the  supersaipts  indicate 
the  phase.  The  chemical  potential  is  usually  written  in  terms  of  the  chemical 
potential  in  an  arbitrary  standard  state,  denoted  by  the  superscript  zero. 


M  -  M®  +  RT  In  p/p® .  (2) 

For  an  ideal  gas  mixture, 

Ml  =  Mi°  +  RT  I®  P/Pi°  .  (3) 


where  pj  is  the  partial  pressure,  equal  to  the  mole  fraction  x|  multiplied  by  P,  the 
total  pressure,  and  the  standard  state  is  usually  pure  component  i. 

For  an  ideal  solid  solution,  the  same  e)q}ression  holds  with  Pi/pi®  replaced  by 
x/xp.  However,  the  standard  state  is  pure  i,  so  xj®  =  1.  The  form  of  eq.  (3)  is  so 
useful  that  it  is  retained  even  for  non-ideal  solutions  with  xj  replaced  by  the 
activity.  a|,  which  may  also  be  considered  a  product  of  x;  multiplied  by  a  non¬ 
ideality  factor,  the  activity  coefficient. 

2. 1  Driving  Force  for  Epitaxy 

For  the  OMVPE  growth  of  GaAs  using  trimethylgallium  (TMGa)  and  arsine 
(ASH3)  the  overall  reaction  is. 

(CH3)3Ga(v)  +  AsH3(v)  =  GaAs(s).  (4) 

Assuming  the  TMGa  and  ASH3  to  completely  decompose  in  the  gas  phase  to 
give  Ga  and  AS4.  an  assumption  to  be  revisited  in  the  discussion  of  kinetics,  to 
follow,  the  reaction  can  be  sinplified: 


V  < 


\H-\ 


if  d  I'or 
uial 


OG 


(5) 


Ga(v)  +  1/4  As4(v)  =  GaAs(s). 

The  equilibrium  condition  is 

p^Ga  +  1/4  pv^s  =  P^GaAs.  (6) 

^  i/ 

M°''Ga  pov^  +  RT  In  p®GaP®/8  ==  M®®GaAs  +  RT  In  aQaAs*  (7) 
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where  the  superscript  'e*  denotes  the  equilibrium  value  of  partial  pressure. 
Thus, 

®GaAs^  P®Ga  P®/^^  ~  ^Ga/^s«  t®) 

where  K  is  the  equilibrium  constant.  This  is  tfie  basic  law  of  mass  acticm. 

When  the  system  is  not  at  equilibrium,  the  tfiermodynamic  driving  force  to 
restore  equilibrium  is 


Ap  =  p''Ga  +  1/4  pv^s^-  pSQaAs. 

(9) 

Ap  =  RT  In  (PGaP^^/P®GaP®/^)* 

(10) 

This  is  the  driving  force  for  epitaxy.  A  situation  Is  Intentionally  created  where 
higher  than  equilibrium  reactant  vapor  pressures  drive  the  system  to  produce 
the  GaAs  solid  desired.  The  maximum  quantity  of  GaAs  solid  that  can  be 
produced  is  simply  the  amount  (the  supersaturation)  that  would  establish 
equilibrium,  and  is  thus  fundamentally  limited  by  thermodynamics  and  the  total 
amount  of  gas  transported  through  the  OMVPE  reactor. 

For  the  OMVPE  growth  of  GaAs  using  arsine  and  TMGa,  the  thermodynamic 
driving  force  at  100  K  is  approximately  80  kcai/mol[1].  This  is  due  to  the 
instability  of  both  arsine  and  TMGa  at  1000  K.  Thermodynamics  also  gives  a 
good  description  of  solid  composition  for  the  OMVPE  growth  of  lll/V  and  I  I/VI 
alloys[2,3]. 

Ordinarily,  in  the  OMVPE  system,  the  growth  rate  is  considerably  less  than  that 
calculated  from  thermodynamics.  Kinetics,  both  surface  reaction  rates  (at  low 
temperatures)  and  diffusion  through  the  gas  phase  (at  higher  temperatures),  are 
not  rapid  enough  to  allow  equilibrium  to  be  established  throughout  the  system 
at  all  times.  This  situation  is  illustrated  by  Fig.  la,  where  Ap  from  eq.  (9)  is 
plotted  versus  reaction  coordinate.  This  allows  the  schematic  represeritation  of 
the  overall,  thermodynamic  driving  force  for  the  growth  reaction,  represented  as 
Ap*.  The  superscript  '* '  denotes  the  chemical  potential  in  the  input  gas  phase, 
where  for  all  reactants  pj  =  p*j.  The  growth  rate  is  proportional  to  the  flux  of 
atoms  diffusing  through  the  boundary  layer,  which  is  identical  to  the  flux  of 
atoms  crossing  the  Interface  Into  the  solid.  The  diagram  shows  schematically 
^e  driving  forces  necessary  to  sustain  this  flux  for  the  diffusion  process  (App) 

and  the  surface  reactions  (/^Js)- 
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Even  in  cases  with  a  large  supersaturation  in  the  input  vapor  phase,  i.e., 
Ap*»0,  near  equilibrium  conditions  may  exist  at  the  growing  solid  surface.  This 
simply  requires  that  the  interface  kinetics  be  much  more  rapid  than  the  diffusion 
kinetics,  i.e.,  the  two  processes  proceed  at  the  same  rate  with  Aps«Ap0.  TNs 
situation,  termed  diffusion  limited  growth,  is  shown  schematically  in  Fig.  1b. 
Using  ordinary  growth  conditions,  with  temperatures  between  approximately 
550  and  800  C,  this  is  the  normal  situation  for  the  OMVPE  growth  of  GaAs,  as 
deduced  from  ^e  nearly  temperature  independent  growth  rate,  as  seen  in  Fig. 
2.  For  surface  kinetically  limited  processes,  tiie  growth  rate  increases 
exponentially  with  Increasing  temperature[4.51.  This  occurs  for  the  OMVPE 
growth  of  GaAs  at  temperatures  below  approximately  550  C  when  TMGa  Is  the 
Ga  precursor,  but  this  temperature  depends  on  the  group  III  precursor  used,  as 
seen  in  Fig.  2  and  discussed  in  Section  3,  below. 

In  the  diffusion  limited  case,  illustrated  schematically  in  Fig.  1b.  the  interfadal 
partial  pressures,  pij,  nearly  satisfy  the  equilibrium  relationship. 

^GaAs/P^GaP'As^'^^aAs-  ('•'•) 

Since  the  input  vapor  is  highly  supersaturated. 

P‘GaP*As>>P'GaP‘Xl:  (12) 

This  is  equivalent  to  stating  that  Ap*»0.  For  the  typical  case 

P*Ga«  4p*As^  (13) 

i.e.,  the  V/lll  ratio  is  »1 ,  as  will  be  discussed  below.  This  means  that  the  Ga  is 
nearly  depleted  at  the  interface. 


P'Ga«P*Ga.  (14) 

while  the  As4  partial  pressure  is  hardly  diminished, 

P'As*P‘As.  (15) 

since  the  same  number  of  As  and  Ga  atoms  are  removed  from  the  vapor  phase 
to  produce  GaAs.  This  situation  makes  the  analysis  of  growth  rate  and  solid 
composition  particularly  simple. 

The  growth  rate  is  proportional  to  the  flux  of  Ga  and  As  atoms  diffusing  through 
the  vapor  to  the  growing  interface.  For  sinr^lidiy,  this  can  be  analyzed  in  terms 
of  diffusion  through  a  boundary  layer  of  thickness  d.  A  more  complete 
description  is  given  in  ttie  following  paper[6].  The  two  fluxes  are  equal,  since 
stoichiometric  GaAs  is  the  only  product.  The  flux  may  be  expressed. 
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(16) 


J  =  DGa(P*Ga  *  p'GaWd. 

where  DQa  'S  the  diffusion  coefficient  of  Ga.  in  whatever  form  it  may  appear 
while  diffusing  through  the  boundary  layer.  In  light  of  eq.  (14),  the  Ga  flux  and 
the  GaAs  growth  rate  are  proportional  to  p*Ga>  observed  6xperimentally[5]. 
Equally  clear  is  that  the  ratio  of  the  concentrations  of  A  and  B  for  alloys  wi^ 
mixing  on  the  group  III  sublattice,  Af.xBxC.  will  be  the  same  as  the  ratio 
p*A/p*B.  assuming  the  diffusion  coefficients  for  the  A  and  B  species  are  nearly 
the  same{7,81. 

2.2  Solution  Thermodynamics 

The  condition  for  thermodynamic  equilibrium  is  expressed  by  eq.  (1)  and  the 
discussion  in  Section  2.1.  Using  these  concepts,  applied  to  the  solid-vapor 
equilibria  of  concern  for  OMVPE,  we  can  calculate  the  composition  of  a 
multicomponent  solid  alloy  from  the  temperatijre  and  the  concentrations  of  the 
various  components  in  the  vapor  phase.  Deviations  from  Ideality  for  the  vapor 
phase  are  commonly  neglected.  However,  non-ideality  in  the  solid  phase  must 
be  considered.  Fortunately,  for  semiconductor  systems  the  solid  can  often  be 
described  using  either  the  regular  solution(9]  or  'delta-lattice-parameter' 
(DLP)[10]  model.  In  both  cases  the  distribution  of  elements  on  a  sublattice  is 
considered  to  be  random,  thus  the  entropy  of  mixing  for  a  pseudobinary  solution 
of  the  type  Ai-xBxC  is  simply  the  ideal  configurational  entropy  of  mixing, 

ASM  =  .R(x  Inx  +  (1-x)  In(l-x)].  (17) 

For  the  regular  solution  model,  the  enthalpy  of  mixing  is  obtained  by  summing 
nearest-neighbor  bond  energies,  yielding. 

AHM  =  x(1-x)QS,  (18) 

where  QS  js  the  interaction  parameter.  The  activity  coefficient  may  be  written, 

lnyj  =  (1-Xi)2Q/RT.  (19) 

Physically,  the  regular  solution  model  cannot  provide  an  accurate,  predictive 
description  of  the  enthalpy  of  mixing  in  semiconductor  alloys.  However,  simple 
models  developed  to  interpret  the  band  gap  and  optical  properties  can  be  used 
to  treat  the  boiling  in  semiconductor  alloys.flO]  The  DLP  model  allows 
accurate  calculation  of  Cfi  in  terms  of  the  difference  in  lattice  parameters 
between  AC  and  BC: 

QS  =  5x10^  (aAC  -  aBC)^(  ^AC  +  aBCV2H  ®.  (20) 

This  first-order  treatment  of  the  enthalpy  of  mixing  is  apparently  equivalent  to 
simply  the  strain  energy  caused  by  tiie  lattice  parameter  difference[1 1].  The 
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solutions  are  predicted  to  be  ideal  for  alloys  from  compounds  with  the  same 
lattice  constant  such  as  GaAs  and  AlAs.  and  to  have  i^sitive  deviations  from 
ideality  for  all  other  alloys.  The  DLP  m(^ei  predicts  a  large  positive  enthalpy  of 
mixing  for  systems  with  a  large  difference  in  lattice  constant.  This  can 
overwhelm  the  negative  configurational  enti-opy  of  mixing  for  temperatures 
below  the  critical  temperature,  Tc,  resulting  in  a  free  energy  versus  corrposition 
curve  with  an  upward  bowing  in  the  center.  This  dictates  ^at  at  equilibrium  a 
random  alloy  In  a  certain  composition  range  will  decompose  into  a  mixture  of 
two  phases,  i.e.,  the  phase  diagram  contains  a  miscibility  gap. 


The  equilibrium  conditions  for  the  ternary(or  pseudobinary)  system  may  be 
obtain^  in  exactly  the  same  way  as  describe  for  binary  systems  in  section  2.1, 
by  equating  the  chemical  potentials  of  the  2  components  in  the  2  phases: 

M''a  +  P''C  =  M®AC.  (21a) 

P''b  +  M''c  =  P^BC-  (21b) 

This  leads  to  two  mass  action  expressions,  similar  to  eq.  (11).  As  discussed  in 
section  2. 1 ,  equilibrium  is  assumed  to  be  established  at  the  interface. 

As  an  example  of  the  use  of  such  calculations  to  understand  epitaxial 
processes,  consider  the  OMVPE  growth  of  GaAs-j.xSbx.  The  2  mass  action 
expressions,  one  for  GaAs  and  one  for  GaSb,  are  solved  simultaneously  with  2 
conservation  equations,  one  for  solid  stoichiometry  and  one  for  solid 
composltion[121.  Complete  pyrolysis  of  the  source  molecules  Is  normally 
assumed.  This  assumption  is  incorrect  for  very  stable  molecules  at  all 
temperatures  and  for  all  molecules  at  very  low  terrperatures,  as  will  be 
discussed  below  in  the  Kinetics  Section.  The  activity  coefficients  of  GaAs  and 
GaSb  in  the  solid  are  calculated  as  described  above  using  the  DLP  model. 

The  calculation  can  be  performed  with  no  adjustable  parameters,  yielding 
solid  composition  versus  vapor  composition  arxi  substrate  temperature  during 
growth.  The  calculated  results  are  compared  with  experimental  data  in  Fig.  3. 
Several  important  aspects  of  VPE  are  illustrated  in  this  rather  complex  figure. 
First,  consider  the  open  data  points,  obtained  for  an  input  V/lil  ratio  (tiie  ratio  of 
the  input  group  V  to  group  III  molar  flow  rates)  of  2.0.  Notice  that  the  calculated 
curve  for  V/lll  =  2.0  fits  the  data  well.  The  Sb  distribution  coefficient,  defined  as 
kSb  =  xSstA''Sb.  whei'e  =  P*TMSb^(P*TMSb+P*AsH  ).  »s  seen  to  be  less 
than  unity.  GaAs  is  more  stable  than  GaSb.  thus  As  is  more  likely  to  bond  to  the 
Ga  on  the  surface  arxi  be  incorporated  into  the  solid.  The  excess  Sb  evaporates 
from  the  surface.  An  additional  implant  point  is  that  the  calculation  for  a  V/lll 
ratio  of  less  ttian  unity  yields  an  antimony  distribution  coefficient  of  unity.  As 
discussed  in  Section  2.1  for  the  case  of  alloys  with  mixing  on  the  group  III 
sublattice,  when  V/lli>1.  essentially  all  of  the  group  III  elements  reaching  the 
interface  are  lrxx)rporated.  The  case  of  GaA^b  with  mixir^  on  the  group  V 
sublattice  with  V/lllsl  is  completely  analogous.  The  establishment  of 
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equilibrium  at  the  interface  while  the  irput  vapor  is  highly  supersaturated 
requires  that  the  group  V  elements  be  virtually  exhausted  at  ttie  interface.  A  final 
point  relative  to  Fig.  3  is  the  presence  of  a  two  solid  phase  region  or  miscibility 
gap.  Because  of  the  large  difference  in  lattice  constant  between  GaAs  and 
GaSb  a  miscibility  gap  exists[13].  However,  when  the  V/lll  ratio  is  less  than 
unity,  the  As  and  Sb  atoms  arriving  in  a  random  pattern  at  the  surface  do  not 
have  time  to  redistribute  themselves  into  GaAs  and  GaSb  rich  areas  before 
being  covered  over  by  the  next  layer.  Thus,  we  are  able  to  grow  metastable 
GaAsi.xSbx  alloys  throughout  the  entire  range  of  solid  composition  as  shown 
by  the  solid  data  points  in  Fig.  3.  Other,  even  less  stable  alloys,  such  as  GaPSb 
and  InPSb  can  also  be  grown  in  this  way[14]. 

2.4  Solid  Phase  Immiscibility 

For  GaAsSb,  tiie  value  of  Tg  Is  approximately  745  C[12].  At  typical  growth 
temperatures,  the  solid  compositions  inside  the  miscibility  gap,  which  covers 
nearly  the  entire  composition  range,  cannot  be  grown  by  liquid  phase  epitaxy 
(LPE)[15].  We  have  already  discussed  the  ability  to  grow  the  metastabte  alloys 
by  OMVPE.  They  can  also  be  grown  by  molecular  beam  epitaxy  (MBE)[16]. 
Recently,  it  has  been  discovered  that  tiiese  alloys  may  also  exhiM  an  ordered, 
superlattice  structure[17],  consisting,  in  the  ideal  case,  of  alternating 
monolayers  of  GaAs  and  GaSb. 

Ordering  in  a  thermodynamic  system  where  the  random  alloy  exhibits  a  large 
positive  enthalpy  of  mixing  is  not  possible  for  a  regular  solution[9],  anotiier 
indication  that  this  model  is  not  physically  appropriate  for  I  HA/  semicorxiuctor 
alloys.  Interestingly,  the  driving  force  for  ordering  is  a  reduction  of  strain  energy, 
the  same  factor  that  leads  to  a  large  positive  enthalpy  of  mixing.  The  ordered 
structures  are  able  to  accommodate  the  two  dissimilar  bond  lengths  when  Aa^O 
with  reduced  distortion  of  the  bonds. 

Ordering  has  now  been  observed  in  essentially  all  lll/V  alloys  grown  by  OMVPE 
and  MBE[1 1  j.  The  {1 1 1}  ordered  structure  (Cu-R)  with  4  variants, 
corresponding  to  the  4  crystallographic  distirx:t  {111}  planes  in  a  cubic  lattice,  is 
normally  observed  for  lll/V  alloys.  Only  2  of  the  variants  are  observed  during 
OMVPE  growth  for  (OOI)-oriented  substrates.  This  is  apparently  due  to  the  lower 
symmetry  of  the  reconstructed,  As-rich  surface.  For  grov^  systems  where 
reconstruction  is  believed  absent,  such  as  LPE,  no  ordering  occurs.  These  arxi 
other  experimental  observations  indicate  that  the  ordering  phenomena  cannot 
be  totally  explained  by  thermodynamic  factors:  Kinetic  factors  appear  to  be 
important. 

3.  KINETICS 

3.1  Survey  of  Precursor  Molecules  Used  for  OMVPE 

Before  discussing  the  chemical  reactions  occurring  during  OMVPE,  it  is 
worthwhile  to  briefly  describe  the  types  of  precursor  molecules  used.  The 
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precursors  typically  contain  the  desired  group  III  or  group  V  element  plus  one  or 
more  of  the  radicals  H.  Cl,  CH3  (methyl),  C2H5  (ethyl),  C3H7  (propyl),  C4H9 
(butyl),  and  others[18].  We  have  recently  entered  the  period  of  ‘designer* 
precursors,  where  special  precursor  molecules  are  designed  for  particular 
applications.  Examples  are  diethylgaltiumchloride  (DEGaCI)  useful  for  selected 
area  growth(19]  and  atomic  layer  epitaxy  (ALE)I20],  tert'arybutylarsine  (TBAs), 
a  less  hazardous  substitute  for  arsine, [21  j  and  trimethylaminealiane  (TMAA), 
an  aluminum  source  that  gives  less  carbOT  contamination  than  the  standard 
trimethylaluminum  (TMAI)[22].  The  groi^)  III  hydrides  would  be  very  desirable 
sources  for  vapor  growth  since  organic  radicals  frequently  result  in  carbon 
contamination  of  the  epitaxial  layers,  as  described  below.  However,  the  group  III 
hydrides  are  too  unstable.  TMAA  was  designed  to  resolve  this  problem  by 
making  an  adduct  compound  with  trimethylamine  to  stabilize  the  alane. 

The  first  requirement  of  a  precursor  is  that  it  be  sufficiently  volatile  to  allow 
acceptable  epitaxial  growth  rates.  The  vapor  pressures,  along  with  other 
properties,  are  tabulated  for  group  V  precursors  in  Table  I.  The  second 
requirement  is  that  heating  to  the  desired  growth  temperature  cause  pyrolysis. 
The  percent  decomposition  versus  temperature  results  obtained  using  an 
isothermal,  flow-tube  reactorI231  for  several  precursors  are  compared  in  terms 
of  the  values  of  T50,  the  temperature  for  50%  pyrolysis  under  the  specific 
conditions  used.  In  Table  I.  The  relative  values  of  T50  are  often  a  reflection  of 
the  strength  of  the  weakest  metal-radical  bond  in  the  precursor  molecuie[18]. 
Arsine  (ASH3)  and  trimethylarsine  (TMAs)  pyrolyze  at  the  highest  temperature 
since  both  the  H-As  and  methyl-As  bonds  are  relatively  strong.  TBAs 
decomposes  at  temperatures  more  than  200°C  less  than  for  arsine  because  of 
the  weak  t-butyl-As  bond.  The  nature  of  the  radicals  is  also  the  most  Important 
factor  determining  carbon  incorporation.  Methyl  radicals  are  quite  reactive,  so 
frequently  lead  to  carbon  incor|X)ration  into  ftie  growing  solid,  especially  for  ^ 
containing  materiais[24].  Less  reactive  radicals,  such  as  t-C4Hg,  apparently 
result  in  very  little  carbon  contamination.  Ethyl  radicals  are  also  much  less  likely 
than  CH3  to  introduce  carbon  into  the  solid[25].  Cl  doping  is  apparently  not  a 
problem  for  lll/V  semiconductors.  H  radicals  actually  act  to  reduce  carbon  by 
supplying  the  H  to  convert  methyl  radicals  to  virtually  inert  methane[18]. 

The  above  discussion  of  pyrolysis  temperature  and  carbon  contamination  is 
somewhat  naive,  sirKe  it  implicitiy  assumes  tiiat  pyrolysis  occurs  by  homolysis 
(or  heterolysis),  i.e.,  by  the  sequential  elimination  of  radicals  from  the  parent 
molecule  until  the  element  is  incorporated  into  the  solid.  For  this  process,  the 
bond  strength  determines  the  pyrolysis  terr^erature  and  the  nature  of  tiie 
resulting  radicals  largely  determines  the  rate  of  carbon  irx;orporation  into  the 
growing  solid.  However,  the  pyrolysis  processes  may  be  considerably  more 
complex.  As  discussed  in  more  detail  in  ftie  next  section,  other  pyrolysis 
processes  may  occur  more  rapidly  than  homolysis  and  produce  pk)  reactive 
radicals,  ^elimination  reactions[18]  occur  without  radical  production,  so  are 
favorable  for  avoiding  carbon  incorporation,  but  do  not  occur  for  methyl  radicals. 
Pyrolysis  of  triemthylarsine  (TMAs)  is  believed  to  occur  via  hydrogenolysis[26]. 

In  this  process,  an  Ha  molecule  interacts  directly  with  the  TMAs  resulting  in  the 
production  of  CH4  plus  (CH3)2AsH.  The  rate  of  this  process  is  not  determined 
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by  the  As-CHs  bond  strength  and  no  CH3  radicals  are  produced.  The  ambient 
can  also  interact  with  radicals  produced  by  homolysis  to  reduce  carbon 
contamination.  The  highly  reactive  methyl  radicals  react  with  an  H2  ambient  to 
produce  methane  and  reactive,  potentially  benefidal,  H  radicals!  18]. 

The  major  factor  driving  the  initial  stages  of  the  search  for  Improved  precursor 
molecules  was  the  demand  for  less  hazardous  precursors.  Arsine  and 
phosphine  are  extremely  dangerous  because  of  their  toxicity  combined  with 
their  high  vapor  pressures  of  many  atmospheres.  Liquid,  organometallic 
sources  are  approximately  100X  less  hazardous  solely  due  to  their  lower  vapor 
pressures,  resulting  in  slower  dispersal  into  8ie  atmos^ere,  which,  in  turn, 
leads  to  lower  concentrations  in  case  of  an  accidental  release  of  the  precursor 
into  the  atmosphere  of  the  laboratory!  18].  Fortunately,  many  organometallic 
molecules  are  also  considerably  less  toxic  than  the  hydrides.  The  outstanding 
example  is  tertiarybutylphosphine  (TBP),  witti  a  toxicity  orders  of  magnitude 
lower  than  that  for  phosphine.  Available  toxicity  data  are  also  included  in  Table 
I. 

The  leading  candidates  for  replacement  of  the  P  and  As  hydrides  are  TBP  and 
either  TBAs  or  ethylarsine  (EAs).  All  of  these  molecules  yield  fairly  stable 
organic  radicals  during  pyrolysis.  They  also  produce  MH2  radicals,  where  M  is 
either  As  or  P,  that  act  to  remove  organic  radicals  from  the  surface.  Thus,  carbon 
contamination  appears  not  to  be  a  problem.  In  fact,  recent  data  on  the  gro>vth  of 
GaAs  and  AIGaAs  by  OMVPE  using  TBAs  indicates  lower  carbon 
concentrations  than  for  growth  using  arsine[27l. 

3.2  Chemical  Reactions 

The  effects  of  chemical  kinetics  are  invariably  observed  during  epitaxial  growth 
when  the  rates  of  the  chemical  reactions  necessary  for  growth  are  slower  than 
the  mass  tran^rt  rates.  Thus,  the  growth  rate  falls  off  at  low  temperatures,  as 
seen  in  Fig.  2  and  discussed  above.  By  comparing  the  temperature 
dependence  of  GaAs  growth  rates  for  several  Ga  precursors,  as  in  Fig.  2,  we 
observe  that  the  pyrolysis  of  the  group  ill  precursor  is  the  rate  determining  step, 
not  surprising  since  the  V/lii  ratio  is  typically  much  greater  than  ur%.  TMGa, 
with  the  highest  radical-Ga  bond  strength,  gives  the  lowest  growth  rate.  The  Ga- 
ethyl  bond  strength  is  lower,  resulting  in  a  higher  grov^  rate  for  triethylgallium 
(TEGa).  The  precursor  with  the  lowest  bond  stren^,  triisobutylgallium  (TIBGa), 
gives  the  highest  growth  rates  at  low  temperatures.  This  is  an  excellent  example 
of  the  interconnection  between  the  details  of  pyrolysis  of  the  precursor  and 
growth  results.  The  practical  consequence  is  that  efficient  low  temperature 
growth  can  only  be  accomplished  by  designir>g  group  III  precursors  with  lower 
bond  strengths. 

Only  recently  have  the  details  of  the  reactions  for  pyrolysis  of  the  group  III  and 
grotqs  V  precursors,  alone  and  in  combination,  been  skidied  in  detail.  The 
chemical  reactions  occum'ng  during  OMVPE  growth  often  involve  complex 
radical  reactions  occurring  both  in  the  gas  phase  and  on  the  surface.  This  is 
perhaps  the  most  complex  and  least  experimentally  accessible  aspect  of  vapor 


phase  growth.  Two  approaches  have  been  adopted  to  study  these  processes: 
1)  The  study  of  chemical  reactions  using  analytical  techniques  such  as  Infrared 
spectroscopy  and  mass  spectroscopy,  either  by  sanpling  ^e  exhaust  or  by  in- 
situ  sampling.  2)  The  study  of  surface  chemical  reactions  using  sophisticated 
surface  analytical  techniques  that  can  normally  be  applied  only  in  UHV 
chambers.  Each  approach  will  be  discussed  briefly. 

Perhaps  the  simplest  approach  to  studying  both  homogeneous  and 
heterogeneous  chemical  reactions  involves  tiie  use  of  a  mass  spectrometer  to 
sample  the  effluent  gas  from  an  ersatz  reactor,  an  isothemnal  flow-tube 
reactor!  18].  This  approach  has  the  advantage  of  yielding  data  from  which  rate 
constants  can  be  determined.  The  disadvantage  is  that  intermediate  spedes 
cannot  be  detected.  Thus,  the  reaction  medianlsms  can  only  be  determined  by 
a  combination  of  sophisticated  techniques  including  the  use  of  isotopic  labeling 
of  the  reactants,  the  comparison  of  reaction  rates  and  products  in  He,  Ha.  atxi 
Da  ambients,  the  addition  and  removal  of  radicals  from  the  system,  and  the  use 
of  various  times,  temperatures,  and  reacta.it  concentrations.  While  it  can  never 
be  claimed  that  the  reaction  mechanism  has  been  unequivocally  determined, 
the  results  and  consequent  pyrolysis  models  may  be  useful  guides  to  the 
selection  of  the  experimental  conditions  to  be  used  for  OMVPE  growth  as  well 
as  the  design  of  the  precursor  molecules,  themselves. 

The  only  example  of  this  approach  considered  here  will  be  a  brief  examination 
of  the  pyrolysis  of  the  novel,  non-hydride  phosphorus  precursor  TBP.  Simply 
determining  the  percent  decomposition  versus  temperature  for  a  fixed  reaction 
time  but  with  variable  input  TBP  concentrations  reveals  that  the  pyrolysis 
process  is  not  unimolecular,  since  the  decomposition  rate  increases  with 
increasing  reactant  concentration[28].  The  products  obtained  versus  reactor 
temperature  for  various  input  concentrations,  are  shown  in  Fig.  4.  They  dearly 
indicate  that  several  reactions  are  occurring  in  parallel,  since  the  product  is 
mainly  C4H8  for  low  Input  TBP  concentrations  and  C4H10  at  higher 
concentrations.  From  the  product  distribution  a  hypothetical  reaction  scheme 
can  be  suggested,  namely  that  C4H8  is  produced  by  the  two  sequential 
unimolecular  reactions  (22)  arxi  (23).  below,  and  that  C4H10  is  (xoduced  by 
radical  attack  of  the  parent  molecule,  removing  an  H  atom,  the  2nd  order 
reaction  (28). 

C4H9PH2  -  C4Hg+PH2  (22) 

C4H9-C4H8+H  (23) 

C4HgPH2  +  C4Hg  C4H9PH  +  C4H10  .  (24) 

The  rate  of  the  second  order  readion.  producing  C4H10.  increases  more 
rapidly  as  the  input  concentration  inaeases  than  the  rates  of  the  first  order 
reactions  (22)  and  (23).  Since  a  number  of  other  models  would  give  similar 
results,  the  hypothetical  model  has  been  tested  by  using[29]  deuterated  TBPd2. 
The  dominant  produd  at  high  concentrations  beeves  C4H9D.  supporting  the 
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proposed  model.  To  further  test  the  proposed  mechanism.  C4H9  radicals  were 
added  to  the  system[30]  by  adding  azo-t-butane  (ATB)  a  compound  that 
pyrolyzes  yielding  2  t-butyl  radicals  and  Inert  N2  at  temperatures  below  those 
required  for  the  onset  of  reaction  (22),  The  temperature  for  the  onset  of  TBP 
pyrolysis  is  dramatically  reduced  and  now  coinddes  with  that  for  ATB.  Clearly 
the  t-butyl  radicals  cause  the  pyrolysis  of  TBP.  Thus,  the  indirect  evidence, 
taken  together,  makes  a  convincing  case  that  TBP  pyrolysis  occurs  mainly  via 
raactions  (22-24).  The  pyrolysis  mechanism  for  TBAs  may  be  similar,  slrx:e  ATB 
has  been  demonstrated  to  have  a  similar  effect  on  the  pyrolysis  ratelSOJ. 

Understanding  the  pyrolysis  processes  for  the  individual  precursors  is  directly 
relevant  to  uriderstanding  and  controlling  the  vapor  phase  growth  techniques. 
As  an  example,  consider  the  growth  of  GaAs  and  GaP  using  TMGa  plus  either 
TMAs  or  TBP.  The  results  of  systematic  studies  of  the  interactions  between 
TMGa  and  TBP  indicate  that  the  PH2  radicals  generated  during  TBP  pyrolysis 
[reaction  (22)1  attack  TMGa  on  the  semiconductor  surface,  enhancing  their 
pyrolysis  rate  by  removing  CH3  radicals[29.31l.  This  Is  also  an  Important 
process  for  removing  methyl  radicals  from  the  system,  resulting  In  a  reduction  in 
carbon  contamination.  Recent  results  for  GaAs  grown  using  TMGa  demonstrate 
TBAs  pyrolysis  probably  also  generates  AsH2  radicals. 

The  chemical  reactions  occurring  specifically  at  the  semiconductor  surface  have 
been  studied  fairly  intensively  by  using  UHV  surface  spectroscopic  probes 
developed  for  the  study  of  catalysis,  such  as  electron  diffraction,  photoemission, 
Auger  electron  spectroscopy,  electron  energy  loss  spectroscopy,  thermally 
programmed  desorption,  arvd  modulated  beam  mass  spectroscopy 
(MBMS)(32,33].  The  elementary  processes  involved  In  UHV  grovrth  techniques, 
chemical  beam  epitaxy(CBE)  for  example,  such  as  the  heterogeneous 
pyrolysis  of  TEG  a,  have  been  the  subject  of  dozens  of  sfr  xJles  within  the  last  few 
years.  Such  studies  have  the  advantage  of  simplicity.  Thus,  they  may  give  more 
direct,  less  ambiguous  information  about  chemical  reactions  than  the  more 
complex  experiments  described  above  for  combined  group  ill  and  group  V 
precursors  in  a  flowing,  atmospheric-pressure  ambient.  Nevertheless,  ^e 
experimental  corxJitions  are  frequently  very  dissimilar  to  those  used  for  growth, 
even  under  UHV  corxlitlons  where  gas  phase  collisions  do  not  occur.  The  times 
involved  are  frequently  very  different  than  those  important  for  growth  processes. 
In  addition,  the  spedes  on  the  surface  may  interact  during  growth.  Thus,  having 
a  specific  surface  coverage  of  a  single  spedes.  although  convenient 
experimentally,  may  reveal  little  about  the  actual  processes  occum’ng  during 
epitaxial  grov\^  where  the  surface  coverage  may  be  much  different.  The  most 
meaning^l  experiments  are  probably  those  involving  the  simultaneous 
preserx^e  of  both  group  III  precursors,  sudi  as  TEGa,  TMGa,  or  trimethylindium 
(TMIn),  and  group  V  molecules,  such  as  As  dimers  and/or  tetramers. 

Using  these  surface  probes,  the  complex  depertdence  of  growth  rate  on 
temperature  for  TEGa  comtx'ned  with  As2  has  been  carefully  studied  and 
modeled.  Martin  and  Whltehouse[32J  have  devised  the  following  scheme  based 
on  their  MBMS  studies.  As  expected,  at  very  low  temperatures.  TEGa  does  not 
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pyrolyze.  At  temperatures  up  to  350®C,  DEGa  Is  formed,  but  desorbs  without 
further  decomposition.  Thus,  the  growth  rate  remains  nearly  zero.  As  the 
temperature  is  increased  above  SSO^’C.  tiie  growth  rate  inaeases.  since  the 
pyrolysis  can  now  proceed  to  elemental  Ga.  However,  a  progressive  deaease 
in  growth  rate  is  observed  as  the  temperature  Is  raised  above  approximately 
440°C.  due  to  an  increasing  rate  of  desorption  of  DEGa.  This  competes  with  the 
growth  process.  Interestingly,  the  basic  scheme  agrees  with  the  earlier 
suggestions  of  Robertson  et  al[34].  based  on  a  simple,  intuitive  model.  Perhaps 
more  Importantly,  these  surface  spectroscopy  experiments  give  information 
useful  in  understanding  more  complex  problems,  such  as  the  rather  complex 
effects  of  temperature  and  group  III  fluxes  on  the  composition  of  alloys,  i.e.. 
GalnAs  and  AIGaAs[32].  This  type  of  information  is  also  expected  to  be 
extremely  valuable  in  the  design  of  new  group  III  precursors  for  CBE. 

3.3.  Motion  of  Surface  Steps 

Surface  analytical  techniques  that  can  be  applied  during  MBE  growth  ( but  not 
for  OMVPE.  because  of  the  electron-absorbing  atmosphere)  have  given 
considerable  information  about  the  atomic  configurations  on  semiconductor 
surfaces  under  various  conditions.  The  unreconstructed  As-rich  (001)  GaAs 
surface  consists  of  As  atoms  forming  only  two  bends  to  underlying  Ga  atoms. 
Tha  two  dangling  bonds  per  surface  atom  give  this  configuration  a  very  high 
energy.  The  energy  is  reduced  dramatically  by  the  formation  of  dimers  between 
adjacent  As  atoms[35l.  These  dimers  run  in  (1 10]  rows  on  the  surface.  An  atom 
being  attached  to  a  kink  moving  along  the  step,  the  process  leading  to  growth, 
thus  sees  two  types  of  sltes[361.  One  involves  the  formation  of  bonds  to  the 
underlying  lattice  atoms  and  the  other  also  Involves  the  formation  of  a  dimer 
bond.  In  addition,  the  location  of  the  dimers  is  correlated  to  the  size  of  the 
underlying  atoms.  These  two  factors  give  rise  to  the  formation  of  2  particular 
variants  of  the  Cu-R  structure  during  growth  on  (OOI)-oriented 
substrates!  1 1 ,36].  The  model  correctly  predicts  the  formation  of  the  same  two 
variants  for  both  alloys  with  mixing  on  the  group  III  sublattice,  such  as  GalnP, 
and  for  those  with  mixing  on  the  Group  V  sublattice,  such  as  GaAsPl36].  The 
model  also  predicts  that  the  direction  of  [1 10]  step  motion  will  determine  which 
of  the  two  variants  will  be  formed.  This  is  dramatically  confirmed  by  recent 
experiments  where  the  direction  of  step  motion  has  been  control!^  by 
patterning  the  surface  with  (1 10]-oriented  groovesI37,38].  The  GalnP  grown  on 
^e  two  sides  of  tiie  groove  is  produced  by  the  motion  of  steps  in  opposite 
directions[37].  This  results  in  each  half  of  the  groove  being  filled  by  a  single 
domain,  consisting  entirely  of  a  material  ordered  with  the  same  variant  of  the 
Cu-R  structure.  Similar,  but  slightly  more  complex,  results  have  been  obtained 
for  GaAsP[38].  This  Is  a  dramatic  demonstration  of  the  control  of  ordering.  These 
are  the  first  maaoscopic  ordered  domains  ever  produced  in  illAf  alloys.  The 
cross  sectional  size  Is  1/2  the  groove  width  by  the  layer  thickness.  The  length  is 
the  dimension  of  the  sutetrate,  i.e.,  the  length  of  the  groove. 

4.  SUMMARY 
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Thermodynamics  acts  as  the  driving  force  for  epitaxial  growth  and  thus  places 
an  upper  limit  on  the  OMVPE  growth  rate.  If  all  of  the  nutrient  phase  in  ^e 
reactor  were  able  to  equilibrate  with  the  solid  epitaxial  layer  being  grown,  the 
system  would  be  operating  at  the  maximum  possible  growth  rate.  This  is 
virtually  never  the  case  b^ause  of  mass  b^ansport  arid  surface  kinetic 
limitations.  For  ‘normal*  growth  rates  and  temperatures  the  chemical  reaction 
kinetics  are  considerably  more  rapid  than  di^sion.  Thus,  most  of  the  vapor 
phase  supersaturation  is  used  to  drive  the  diffusion  process.  As  a  result,  the 
vapor  phase  adjacent  to  tfie  growing  interface  is  nearly  in  equilibrium  with  the 
solid.  This  situation  allows  the  construction  of  a  simple  model  that  can  be  used 
to  accurately  predict  the  solid  compositions  in  lll/V  alloys  in  terms  of  the  input 
vapor  pressures  and  sutetrate  temperature  during  growth. 

Miscibility  gaps  exist  in  many  lil/V  alloys  due  to  the  large  difference  in  lattice 
parameters  between  the  binary  end  components.  OMVPE  has  been  used  for 
the  growth  of  metastable  alloys.  The  key  to  growing  these  metastable  alloys  is 
the  input  V/lll  ratio.  When  V/lll<1,  all  of  tiie  As  and  Sb  reaching  the  interface  are 
incorporated.  The  random  As  and  Sb  distribution  expected  from  the  random 
impingement  from  the  vapa  phase  might  be  expected  to  result  in  a  nearly 
random  alloy.  However,  ordered  structures  are  commonly  observed  in  electron 
diffraction  patterns.  The  ordering  is  due  to  the  high  strain  energy  inherent  in  a 
random  solid  alloy  consisting  of  bonds  witii  different  preferred  lengths.  The 
formation  of  ordered  structures  reduces  this  miaoscoptc  strain  energy. 

The  kinetic  aspects  of  epitaxy  are  by  far  the  most  complex  and,  consequently, 
the  least  understood.  We  are  beginning  to  understand  the  reactions  involved  in 
the  homogeneous  and  heterogeneous  pyrolysis  of  the  wide  variety  of  group  III 
and  group  V  precursors  used  for  OMVPE.  These  reactions  are  often  complex, 
involving  second  order  processes  such  as  the  attack  of  the  parent  molecules  by 
free  radicals  produced  during  homolysis.  The  processes  occum'ng  on  the 
surface  that  lead  to  growth  are  even  more  complex.  Both  simple  bond  breaking 
to  produce  radicals  and  radical  attack  reactions  apparently  occur 
heterogeneously.  In  addition,  desorption  of  parent  molecules  and  intermediates 
controls  the  growth  rate  under  certain  conditions,  as  convincingly  demonstrated 
for  the  CBE  growth  of  GaAs  from  TEGa  and  cracked  arsine.  The  physical 
processes,  such  as  surface  reconstruction  and  step  motion,  also  have  profound 
effects  on  the  solid  formed,  including  the  occurrence  of  phase  separation  and 
ordering. 

ACKNOWLEDGMENTS 

The  author  gratefully  acknowledges  the  financial  support  of  the  Department  of 
Energy,  the  Office  of  Naval  Research,  and  the  Airforce  Office  of  Scientific 
Research. 


13 


Table  I:  Properties  of  group  V  sources  used  for  epitaxy.  The  values  of  vapor 
pressure  and  toxicity  were  obtained  from  ref.  18,  unless  otherwise  indicated. 
The  values  of  pyrolysis  temperature  (Tso)  were  obtained  from  similar 
experiments  in  an  atmospheric-pressure,  flow-tube  reactor  with  a  residence 
time  of  a  few  seconds. 


Precursor 

Vapor  Pressure 
p.Torr/T.°C 

T5Q(!C1 

Toxicity 

LCSQ"  (TLVl 

Carbon  Doping 
TMGaA-Mln 

PH3 

850 

11-50  (0.3) 

Very  Low 

TBP 

141/10 

450 

>1100 

Low 

AsH3 

600 

5-50  (0.05) 

Very  Low 

TMAs 

238/20 

530 

20,000 

Very  High 

TEAS 

5/20 

490 

500-1000*** 

High 

DMAs 

176/0 

460 

130 

High 

DBAs 

40/20 

440 

300 

Low 

EAs 

145/-7a 

440 

Low 

TBAs 

96/10 

380 

70 

Very  Low 

PhAsH2 

(C6H5)AsH2 

1.8/20^ 

High 

*  TLV  based  on  average  eight  hour  work  day 
**  Based  on  study  of  rat  mortality  after  4  hour  exposure 
***Orai  dose  in  mg  of  material  per  kg  of  animal  weight. 

a)  D.M.  Speckman  and  J.P.  Wendt,  J.  Crystal  Growth  105, 275  (1990). 

b)  A.  Brauers,  O.  Kayser,  R.  Kail,  H.  Heinecke,  P.  Balk,  and  H.  Hofmann,  J. 
Crystal  Growth  93,  7  (1988). 
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FIGURE  CAPTIONS 


Figure  1.  Schematic  diagram  of  chemical  potential  versus  reaction  coordinate, 
showing  the  drop  in  chemical  potential  required  for  each  step  in  the  growth 
sequence  to  keep  all  rates  equal.  The  differences  in  individual  chemical 
potentials  can  alternatively  be  thought  of  as  ratios  of  partial  pressures  of  the 
reactants;  a)  the  general  case  (After  Stringfellowr[39])  and  b)  the  case  of  rapid 
surface  kinetics,  i.e.,  with  (After  Strin^ellow  [40]). 

Figure  2.  Growth  efficiency  (growth  rate/group  III  molar  flow  rate)  versus 
reciprocal  temperature  for  various  Ga  alkyls;  TMGa  W.  TEGa  (•),  and  TIBGa 
(A).  (After  Stringfellow[41]).  The  data  are  from  Plass  et  al[42]. 

Figure  3.  Solid  versus  vapor  composition  for  the  alloy  GaAs|.xSbx.  The  data 
are  from  the  work  of  Cooper  et  al[43]  for  V/lll  =  2.0(0),  and  V/lll  =  0.5((i),  arxl  the 
work  of  Stringfellow  and  Cherng[44i(e).  The  curves  were  calculated  for  various 
V/lll  ratios.  The  broken  sections  of  each  curve  represent  the  calculated  regions 
of  solid  immiscibility.  (After  Stringfellow  and  Cherng[44]). 

Figure  4.  Products  of  TBP  decomposition  versus  temperature  for  several  initial 
TBP  concentrations.  (After  Li  et  al(45]). 
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